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Preface 



This book results from my teaching and research activities at the Technical Uni- 
versity of Hamburg-Harburg (TUHH), Germany. It is based on my German book 
“Mobilfunkkanale — Modellierung, Analyse und Simulation” published by Vieweg & 
Sohn, Braunschweig/Wiesbaden, Germany, in 1999. The German version served as a 
text for the lecture Modern Methods for Modelling of Networks , which I gave at the 
TUHH from 1996 to 2000 for students in electrical engineering at masters level. 

The book mainly is addressed to engineers, computer scientists, and physicists, who 
work in the industry or in research institutes in the wireless communications field and 
therefore have a professional interest in subjects dealing with mobile fading channels. 
In addition to that, it is also suitable for scientists working on present problems of 
stochastic and deterministic channel modelling. Last, but not least, this book also is 
addressed to master students of electrical engineering who are specialising in mobile 
radio communications. 

In order to be able to study this book, basic knowledge of probability theory and 
system theory is required, with which students at masters level are in general 
familiar. In order to simplify comprehension, the fundamental mathematical tools, 
which are relevant for the objectives of this book, are recapitulated at the beginning. 
Starting from this basic knowledge, nearly all statements made in this book are 
derived in detail, so that a high grade of mathematical unity is achieved. Thanks 
to sufficient advice and help, it is guaranteed that the interested reader can verify 
the results with reasonable effort. Longer derivations interrupting the flow of the 
content are found in the Appendices. There, the reader can also find a selection 
of MATLAB-programs, which should give practical help in the application of the 
methods described in the book. To illustrate the results, a large number of figures 
have been included, whose meanings are explained in the text. Use of abbreviations 
has generally been avoided, which in my experience simplifies the readability consid- 
erably. Furthermore, a large number of references is provided, so that the reader is led 
to further sources of the almost inexhaustible topic of mobile fading channel modelling. 

My aim was to introduce the reader to the fundamentals of modelling, analysis, 
and simulation of mobile fading channels. One of the main focuses of this book is 
the treatment of deterministic processes. They form the basis for the development 
of efficient channel simulators. For the design of deterministic processes with given 
correlation properties, nearly all the methods known in the literature up to now 
are introduced, analysed, and assessed on their performance in this book. Further 
focus is put on the derivation and analysis of stochastic channel models as well 
as on the development of highly precise channel simulators for various classes of 
frequency-selective and frequency-nonselective mobile radio channels. Moreover, a 
primary topic is the fitting of the statistical properties of the designed channel models 
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to the statistics of real-worlcl channels. 

At this point, I would like to thank those people, without whose help this book 
would never have been published in its present form. First, I would like to express 
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INTRODUCTION 



1.1 THE EVOLUTION OF MOBILE RADIO SYSTEMS 

For several years, the mobile communications sector has definitely been the fastest- 
growing market segment in telecommunications. Experts agree that today we are 
just at the beginning of a global development, which will increase considerably 
during the next years. Trying to find the factors responsible for this development, 
one immediately discovers a broad range of reasons. Certainly, the liberalization 
of the telecommunication services, the opening and deregulation of the European 
markets, the topping of frequency ranges around and over 1 GHz, improved modu- 
lation and coding techniques, as well as impressive progress in the semiconductor 
technology (e.g., large-scale integrated CMOS- and GaAs-technology) , and, last but 
not least, a better knowledge of the propagation processes of electromagnetic waves 
in an extraordinary complex environment have made their contribution to this success. 

The beginning of this turbulent development now can be traced to more than 40 
years ago. The first generation mobile radio systems developed at that time were 
entirely based on analog technique. They were strictly limited in their capacity of 
subscribers and their accessibility. The first mobile radio network in Germany was in 
service between 1958 and 1977. It was randomly named A-net and was still based on 
manual switching. Direct dialling was at first possible with the B-net, introduced in 
1972. Nevertheless, the calling party had to know where the called party was located 
and, moreover, the capacity limit of 27 000 subscribers was reached fairly quickly. The 
B-net was taken out of service on the 31st of December 1994. Automatic localization 
of the mobile subscriber and passing on to the next cell was at first possible with the 
cellular C-net introduced in 1986. It operates at a frequency range of 450 MHz and 
has a Germany-wide accessibility with a capacity of 750 000 subscribers. 

Second generation mobile radio systems are characterized by digitalization of the 
networks. The GSM standard (GSM: Groupe Special Mobile) 1 developed in Europe 
is generally accepted as the most elaborated standard worldwide. The D-net, brought 
into service in 1992, is based on the GSM standard. It operates at a frequency range 
of 900 MHz and offers all subscribers a Europe-wide coverage. In addition to this, 
the E-net (Digital Cellular System, DCS 1800) has been running parallel to the 



1 By now GSM stands for “Global System for Mobile Communications” . 
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D-net since 1994, operating at a frequency range of 1800 MHz. Mainly, these two 
networks only differ in their respective frequency range. In Great Britain, however, 
the DCS 1800 is known as PCN (Personal Communications Network). Estimates say 
the amount of subscribers using mobile telephones will in Europe alone grow from 92 
million at present to 215 million at the end of 2005. In consequence, it is expected 
that in Europe the number of employees in this branch will grow from 115 000 at 
present to 1.89 million (source of information: Lehman Brothers Telecom Research 
estimates). The originally European GSM standard has in the meantime become 
a worldwide mobile communication standard that has been accepted by 129 (110) 
countries at the end of 1998 (1997). The network operators altogether ran 256 GSM 
networks with over 70.3 million subscribers at the end of 1997 worldwide. But only 
one year later (at the end of 1998), the amount of GSM networks had increased to 
324 with 135 million subscribers. In addition to the GSM standard, a new standard 
for cordless telephones, the DECT standard (DECT: Digital European Cordless 
Telephone), was introduced by the European Telecommunications Standard Institute 
(ETSI) . The DECT standard allows subscribers moving at a fair pace to use cordless 
telephones at a maximum range of about 300 m. 

In Europe, third generation mobile radio systems is expected to be practically ready 
for use at the beginning of the twenty-first century with the introduction of the 
Universal Mobile Telecommunications System (UMTS) and the Mobile Broadband 
System (MBS). With UMTS, in Europe one is aiming at integrating the various 
services offered by second generation mobile radio systems into one universal system 
[Nie92]. An individual subscriber can then be called at any time, from any place 
(car, train, aircraft, etc.) and will be able to use all services via a universal terminal. 
With the same aim, the system IMT 2000 (International Mobile Telecommunications 
2000) 2 is being worked on worldwide. Apart from that, UMTS/IMT 2000 will also 
provide multimedia services and other broadband services with maximum data rates 
up to 2 Mbit/s at a frequency range of 2 GHz. MBS plans mobile broadband services 
up to a data rate of 155 Mbit/s at a frequency range between 60 and 70 GHz. This 
concept is aimed to cover the whole area with mobile terminals, from fixed optical 
fibre networks over optical fibre connected base stations to the indoor area. For 
UMTS/IMT 2000 as for MBS, communication by satellites will be of vital importance. 

From future satellite communication it will be expected besides supplying areas 
with weak infrastructure — that mobile communication systems can be realized 
for global usage. The present INMARSAT-M system, based on four geostationary 
satellites (35 786km altitude), will at the turn of the century be replaced by satellites 
flying on non-geostationary orbits at medium height (Medium Earth Orbit, MEO) 
and at low height (Low Earth Orbit, LEO). The MEO satellite system is represented 
by ICO with 12 satellites circling at an altitude of 10 354 km, and typical represen- 
tatives of the LEO satellite systems are IRIDIUM (66 satellites, 780km altitude), 
GLOBALSTAR (48 satellites, 1414km altitude), and TELEDESIC (288 satellites, 



2 IMT 2000 was formally known as FPLMTS (Future Public Land Mobile Telecommunications 
System) . 
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1400 km altitude) 3 [Pad95]. A coverage area at 1.6 GHz is intended for hand-portable 
terminals that have about the same size and weight as GSM mobile telephones today. 
On the 1st of November 1998, IRIDIUM took the first satellite telephone network 
into service. An Iridium satellite telephone cost 5 999 DM in April 1999, and the price 
for a call was, depending on the location, settled at 5 to 20 DM per minute. Despite 
the high prices for equipment and calls, it is estimated that in the next ten years 
about 60 million customers worldwide will buy satellite telephones. 

At the end of this technical evolution from today’s point of view is the development 
of the fourth generation mobile radio systems. The aim of this is integration of 
broadband mobile services, which will make it necessary to extend the mobile 
communication to frequency ranges up to 100 GHz. 

Before the introduction of each newly developed mobile communication systems a large 
number of theoretical and experimental investigations have to be made. These help to 
answer open questions, e.g., how existing resources (energy, frequency range, labour, 
ground, capital) can be used economically with a growing number of subscribers and 
how reliable, secure data transmission can be provided for the user as cheaply and as 
simple to handle as possible. Also included are estimates of environmental and health 
risks that almost inevitably exist when mass-market technologies are introduced and 
that are only to a certain extent tolerated by a public becoming more and more critical. 
Another boundary condition growing in importance during the development of new 
transmission techniques is often the demand for compatibility with existing systems. 
To solve the technical problems related to these boundary conditions, it is necessary 
to have a firm knowledge of the specific characteristics of the mobile radio channel. 
The term mobile radio channel in this context is the physical medium that is used 
to send the signal from the transmitter to the receiver [Pro95]. However, when the 
channel is modelled, the characteristics of the transmitting and the receiving antenna 
are in general included in the channel model. The basic characteristics of mobile radio 
channels are explained later. The thermal noise is not taken into consideration in the 
following and has to be added separately to the output signal of the mobile radio 
channel, if necessary. 

1.2 BASIC KNOWLEDGE OF MOBILE RADIO CHANNELS 

In mobile radio communications, the emitted electromagnetic waves often do not 
reach the receiving antenna directly due to obstacles blocking the line-of-sight path. In 
fact, the received waves are a superposition of waves coming from all directions due to 
reflection, diffraction, and scattering caused by buildings, trees, and other obstacles. 
This effect is known as multipath propagation. A typical scenario for the terrestrial 
mobile radio channel is shown in Figure 1.1. Due to the multipath propagation, the 
received signal consists of an infinite sum of attenuated, delayed, and phase-shifted 
replicas of the transmitted signal, each influencing each other. Depending on the 
phase of each partial wave, the superposition can be constructive or destructive. Apart 



3 Originally TELEDESIC planned to operate 924 satellites circling at an altitude between 695 and 
705 km. 
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from that, when transmitting digital signals, the form of the transmitted impulse 
can be distorted during transmission and often several individually distinguishable 
impulses occur at the receiver due to multipath propagation. This effect is called the 
impulse dispersion. The value of the impulse dispersion depends on the propagation 
delay differences and the amplitude relations of the partial waves. We will see 
later on that multipath propagation in a frequency domain expresses itself in the 
non-ideal frequency response of the transfer function of the mobile radio channel. 
As a consequence, the channel distorts the frequency response characteristic of the 
transmitted signal. The distortions caused by multipath propagation are linear and 
have to be compensated for on the receiver side, for example, by an equalizer. 




Figure 1.1: Typical mobile radio scenario illustrating multipath propagation in a 

terrestrial mobile radio environment. 



Besides the multipath propagation, also the Doppler effect has a negative influence 
on the transmission characteristics of the mobile radio channel. Due to the movement 
of the mobile unit, the Doppler effect causes a frequency shift of each of the partial 
waves. The angle of arrival a n , which is defined by the direction of arrival of the nth 
incident wave and the direction of motion of the mobile unit as shown in Figure 1.2, 
determines the Doppler frequency ( frequency shift) of the nth. incident wave according 
to the relation 



fn ■ — fmax COS Qt n . 



( 1 . 1 ) 
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In this case, fmax is the maximum Doppler frequency related to the speed of the mobile 
unit v, the speed of light Co, and the carrier frequency fo by the equation 

fmax = - fo ■ (1-2) 

Co 

The maximum (minimum) Doppler frequency, i.e., /„ = fmax ( f n = - fmax), is 
reached for a n = 0 (a n = i r). In comparison, though, /„ = 0 for a n = 7t/2 and 
a n = 37t/2. Due to the Doppler effect, the spectrum of the transmitted signal 
undergoes a frequency expansion during transmission. This effect is called the 
frequency dispersion. The value of the frequency dispersion mainly depends on the 
maximum Doppler frequency and the amplitudes of the received partial waves. In the 
time domain, the Doppler effect implicates that the impulse response of the channel 
becomes time- variant. One can easily show that mobile radio channels fulfil the 
principle of superposition [Opp75, Lue90] and therefore are linear systems. Due to 
the time- variant behaviour of the impulse response, mobile radio channels therefore 
generally belong to the class of linear time- variant systems. 




Figure 1.2: 



Angle of arrival a„ of the nth incident wave illustrating the Doppler effect. 



Multipath propagation in connection with the movement of the receiver and/or the 
transmitter leads to drastic and random fluctuations of the received signal. Fades of 
30 to 40 dB and more below the mean value of the received signal level can occur 
several times per second, depending on the speed of the mobile unit and the carrier 
frequency [Jak93]. A typical example of the behaviour of the received signal in mobile 
communications is shown in Figure 1.3. In this case, the speed of the mobile unit 
is v = 110 km/h and the carrier frequency is / 0 = 900 MHz. According to (1.2), 
this corresponds to a maximum Doppler frequency of fmax = 91 Hz. In the present 
example, the distance covered by the mobile unit during the chosen period of time 
from 0 to 0.327s is equal to 10 m. 

In digital data transmission, the momentary fading of the received signal causes burst 
errors , i.e., errors with strong statistical connections to each other [Bla84] . Therefore, 
a fading interval produces burst errors, where the burst length is determined by the 
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Time, t (s) 



Figure 1.3: Typical behaviour of the received signal in mobile communications. 



duration of the fading interval for which the expression duration of fades has been 
introduced in [Kuc82]. Corresponding to this, a connecting interval produces a bit 
sequence almost free of errors. Its length depends on the duration of the connecting 
interval for which the term connecting time interval has been established [Kuc82]. 
As suitable measures for error protection, high performance procedures for channel 
coding are called in to help. Developing and dimensioning of codes require knowledge 
of the statistical distribution of the duration of fades and of the connecting time 
intervals as exact as possible. The task of channel modelling now is to record and to 
model the main influences on signal transmission to create a basis for the development 
of transmission systems [Kit82] . 

Modern methods of modelling mobile radio channels are especially useful, for they 
not only can model the statistical properties of real-world (measured) channels 
regarding the probability density function (first order statistics) of the channel 
amplitude sufficiently enough, but also regarding the level-crossing rate (second order 
statistics) and the average duration of fades (second order statistics). Questions 
connected to this theme will be treated in detail in this book. Mainly, two goals 
are aimed at. The first one is to find stochastic processes especially suitable for 
modelling frequency-nonselective and frequency-selective mobile radio channels. In 
this context, we will establish a channel model described by ideal (not realizable) 
stochastic processes as the reference model or as the analytical model. The second 
goal is the derivation of efficient and flexible simulation models for various typical 
mobile radio scenarios. Following these aims, the relations shown in Figure 1.4, which 
demonstrates the connections between the physical channel, the stochastic reference 
model, and the therefrom derivable deterministic simulation model, will accompany 
us throughout the book. The usefulness and the quality of a reference model and 
the corresponding simulation model are ultimately judged on how well its individual 
statistics can be adapted to the statistical properties of measured or specified channels. 
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Figure 1.4: Relationships between the physical channel, the stochastic reference model, 

the deterministic simulation model, and the measurement or specification. 



1.3 STRUCTURE OF THIS BOOK 

A good knowledge of statistics and system theory are the necessary tools for engineers 
in practice as well as for scientists working in research areas, making the approach to 
a deeper understanding of channel modelling possible. Therefore, in Chapter 2 some 
important terms, definitions, and formulae often referred to in following chapters will 
be recapitulated. Chapter 2 makes the reader familiar with the nomenclature used in 
the book. 

Building on the terms introduced in Chapter 2, in Chapter 3 Rayleigh and Rice 
processes are dealt with as reference models to describe frequency-nonselective 
mobile radio channels. These are at first described in general (Section 3.1). Then, 
a description of the most frequently used Doppler power spectral densities (Jakes 
or Clarke power spectral density and Gaussian power spectral density) and their 
characteristic quantities such as the Doppler shift and the Doppler spread are given 
(Section 3.2). After that, in Section 3.3 the statistical properties of the first kind 
(probability density of amplitude and phase) and of the second kind (level-crossing 
rate and average duration of fades) are investigated. Chapter 3 ends with an analysis 
of the statistics of the duration of fades of Rayleigh processes. 

In Chapter 4, it is at first made clear that, from the developers of simulation models 
point of view, analytical models represent reference models to a certain extent. 
Their relevant statistical properties will be modelled sufficiently exactly with the 
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smallest possible realization expenditure. To solve this problem, various deterministic 
and statistic methods have been proposed in the literature. The heart of many 
procedures for channel modelling is based on the principle that filtered Gaussian 
random processes can be approximated by a sum of weighted harmonic functions. 
This principle in itself is not at all new, but can historically be traced back to 
basic works of S. O. Rice [Ric44, Ric45]. In principle, all attempts to compute the 
parameters of a simulation model can be classified as either statistic, deterministic 
or as a combination of both. A fact though is that the resulting simulation model 
is definitely of pure deterministic nature, which is made clear in Section 4.1. The 
analysis of the elementary properties of deterministic simulation systems is therefore 
mainly performed by the system theory and signal theory (Section 4.2). Investigating 
the statistical properties of the first kind and of the second kind, however, we will 
again make use of the probability theory and statistics (Section 4.3). 

Chapter 5 contains a comprehensive description of the most important procedures 
presently known for computing the model parameters of deterministic simulation 
models (Sections 5.1 and 5.2). The performance of each procedure will be assessed 
with the help of quality criteria. Often, also the individual methods are compared in 
their performance to allow the advantages and disadvantages stand out. Chapter 5 
ends with an analysis of the duration of fades of deterministic Rayleigh processes 
(Section 5.3). 

It is well known that the statistics (of the first kind and of the second kind) of 
Rayleigh and Rice processes can only be influenced by a small number of parameters. 
On the one hand, this makes the mathematical description of the model much easier, 
but on the other hand, however, it narrows the flexibility of these stochastic processes. 
A consequence of this is that the statistical properties of real-world channels can 
only be roughly modelled with Rayleigh and Rice processes. For a finer adaptation 
to reality, one therefore needs more sophisticated model processes. Chapter 6 
deals with the description of stochastic and deterministic processes for modelling 
frequency-nonselective mobile radio channels. The so-called extended Suzuki processes 
of Type I (Section 6.1) and of Type II (Section 6.2) as well as generalized Rice and 
Suzuki processes (Section 6.3) are derived and their statistical properties are analysed. 
Apart from that, in Section 6.4, a modified version of the Loo model is introduced, 
containing the classical Loo model as a special case. To demonstrate the usefulness 
of all channel models suggested in this chapter, the statistical properties (probability 
density of the channel amplitude, level-crossing rate and average duration of fades) 
of each model are fitted to measurement results in the literature and are compared 
with the corresponding simulation results. 

Chapter 7 is dedicated to the description of frequency-selective stochastic and 
deterministic channel models beginning with the ellipses model introduced by Parsons 
and Bajwa, illustrating the path geometry for multipath fading channels (Section 7.1). 
In Section 7.2, a description of linear time-variant systems is given. With the help 
of system theory, four important system functions are introduced allowing us to 
describe the input-output behaviour of linear time-variant systems in different 
ways. Section 7.3 is devoted to the theory of linear time- variant stochastic systems 
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going back to Bello [Bel63]. In connection with this, stochastic system functions 
and characteristic quantities derivable from these are defined. Also the reference 
to frequency-selective stochastic channel models is established and, moreover, the 
channel models for typical propagation areas specified in the European work group 
COST 207 [COS89] are given. Section 7.4 deals with the derivation and analysis 
of frequency-selective deterministic channel models. Chapter 7 ends with the de- 
sign of deterministic simulation models for the channel models according to COST 207. 

Chapter 8 deals with the derivation, analysis, and realization of fast channel 
simulators. For the derivation of fast channel simulators, the periodicity of harmonic 
functions is exploited. It is shown how alternative structures for the simulation of 
deterministic processes can be derived. In particular, for complex Gaussian random 
processes it is extraordinarily easy to derive simulation models merely based on 
adders, storage elements, and simple address generators. During the actual simulation 
of the complex- valued channel amplitude, time-consuming trigonometric operations as 
well as multiplications are then no longer required. This results in high-speed channel 
simulators, which are suitable for all frequency-selective and frequency-nonselective 
channel models dealt with in previous chapters. Since the proposed principle can be 
generalized easily, we will in Chapter 8 restrict our attention to the derivation of fast 
channel simulators for Rayleigh channels. Therefore, we will exclusively employ the 
discrete-time representation and will introduce so-called discrete-time deterministic 
processes in Section 8.1. With these processes there are new possibilities for indirect 
realization. The three most important of them are introduced in Section 8.2. In the 
following Section 8.3, the elementary and statistical properties of discrete deter- 
ministic processes are examined. Section 8.4 deals with the analysis of the required 
realization expenditure and with the measurement of the simulation speed of fast 
channel simulators. Chapter 8 ends with a comparison between the Rice method and 
the filter method (Section 8.5). 




2 



RANDOM VARIABLES, 
STOCHASTIC PROCESSES, 
AND DETERMINISTIC 
SIGNALS 



Besides clarifying the used nomenclature, we will in this chapter introduce some 
important terms, which will later often be used in the context of describing stochastic 
and deterministic channel models. However, the primary aim is to familiarize the 
reader with some basic principles and definitions of probability, random signals, 
and systems theory, as far as it is necessary for the understanding of this book. 
A complete and detailed description of these subjects will not be presented here; 
instead, some relevant technical literature will be recommended for further studies. 
As technical literature for the subject of probability theory, random variables, and 
stochastic processes, the books by Papoulis [Pap91], Peebles [Pee93], Therrien [The92], 
Dupraz [Dup86], as well as Shanmugan and Breipohl [Sha88] are recommended. 
Also the classical works of Middleton [Mid60], Davenport [Dav70], and the book 
by Davenport and Root [Dav58] are even nowadays still worth reading. A modern 
German introduction to the basic principles of probability and stochastic processes 
can be found in [Boe98, Bei97, Hae97]. Finally, the excellent textbooks by Oppenheim 
and Schafer [Opp75], Papoulis [Pap77], Rabiner and Gold [Rab75], Kailath [Kai80], 
Unbehauen [Unb90], SchiiBler [Sch91], and Fettweis [Fet96] provide a deep insight into 
systems theory as well as into the principles of digital signal processing. 

2.1 RANDOM VARIABLES 

In the context of this book, random variables are of central importance, not only to the 
statistical but also to the deterministic modelling of mobile radio channels. Therefore, 
we will at first review some basic definitions and terms which are frequently used in 
connection with random variables. 

An experiment whose outcome is not known in advance is called a random experiment. 
We will call points representing the outcomes of a random experiment sample points 
s. A collection of possible outcomes of a random experiment is an event A. The 
event A = {s} consisting of a single element s is an elementary event. The set of 
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all possible outcomes of a given random experiment is called the sample space Q of 
that experiment. Hence, a sample point is an element of the event, i.e., s G A, and the 
event itself is a subset of the sample space, i.e., A c Q. The sample space Q is called 
the certain event , and the empty set or null set, denoted by 0, is the impossible event. 
Let A be a class (collection) of subsets of a sample space Q. In probability theory, A 
is often called a-field (or a-algebra), if and only if the following conditions are fulfilled: 



(i) The empty set 0 G A. 

(ii) If A G A, then also Q — A G A, i.e., if the event A is an element of the class A , 
then so is its complement. 

(iv) If A n G A (n = 1,2, . . .), then also U ^LiA n G A , i.e., if the events A n are all 
elements of the class A, then so is their countable union. 

A pair ( Q,A ) consisting of a sample space Q and a er-field A is called a measurable 
space. 

A mapping P : A — > IR is called the probability measure or briefly probability, if the 
following conditions are fulfilled: 

(i) If A € A, then 0 < P(A) < 1. 

(ii) P(Q) = 1. 

(iii) If A n G A (n = 1,2,...) with U^^An G A and A n fl Ak = 0 for any n^k, then 
also P(U~ =1 A„) = £~iP(A„). 

A probability space is the triple (Q, A,P). 

A random variable p £ Q is a mapping which assigns to every outcome s of a random 
experiment a number p(s), i.e., 

p : Q — > IR , s i — > p(s) . (2-1) 

This mapping has the property that the set {s|/z(s) < x} is an event of the considered 
cr-algebra for all x G IR, i.e., {s|/i(s) < i} G A Hence, a random variable is a function 
of the elements of a sample space Q. 

For the probability that the random variable p is less or equal to x, we use the 
simplified notation 

P{p < x) := P({s\p(s) < x}) (2-2) 

in the sequel. 

Cumulative distribution function: The function F ;i , defined by 

: IR. -> [0, 1] , xi — > F^{x) = P{p < x ) , (2.3) 

is called the cumulative distribution function of the random variable p. The cumulative 
distribution function F^{x) satisfies the following properties: a) P p (— oo) = 0; b) 
F)j,( oo ) = 1; and c) F^(x) is non-decreasing, i.e., P^aq) < P^aq) if X\ < aq. 
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Probability density function: The function p^, defined by 

Pn : IR . 1R , x i — » Pn(x) = ’ ( 2 -4) 

is called the probability density function (or probability density or simply density) of 
the random variable /i, where it is assumed that the cumulative distribution function 
Ffj(x) is differentiable with respect to x. The probability density function p^{x) 
satisfies the following properties: a) p^x) > 0 for all x; b) fT^Pni 3 ') dx = ^ and c) 

F n(x) = f-ooPn( x ) dx - 

Joint cumulative distribution function: The function F' /i|/j2 , defined by 

F MlM2 : IR 2 -> [0, 1] , (xi,x 2 )> — » F flltl2 (xi,X 2 ) = P(pi < xi,p 2 < x 2 ) , (2.5) 

is called the joint cumulative distribution function (or bivariate cumulative distribution 
function ) of the random variables p\ and p 2 . 

Joint probability density function: The function p MlM2 , defined by 

Pmn 2 '■ ^ JR. {xi,x 2 )* — > p f _ lll _ l2 (x 1 ,x 2 ) = 9 ’ ( 2 - 6 ) 

is called the joint probability density function (or bivariate density function or simply 
bivariate density) of the random variables p\ and p 2 , where it is assumed that the joint 
cumulative distribution function F Ml/l2 ( x±,x 2 ) is partially differentiable with respect 
to X\ and x 2 . 

The random variables p\ and p 2 are said to be statistically independent , if the events 
{s|/xi(s) < Xi} and {s|/i 2 (s) < x 2 } are independent for all # 1 , 2:2 £ IR. In this case, 
we can write -F Ml|U2 (a; i,x 2 ) = F^(x i) • F t _ l2 (x 2 ) and p Ml([t2 (a: i,x 2 ) = p^{x i) -p M2 (x 2 ). 

The marginal probability density functions (or marginal densities) of the joint 
probability density function p l i lll2 (xi,x 2 ) are obtained by 



OO 



Pm (*l) = 


J PuiV2( x h X 2) dx 2 , 

— 00 


(2.7a) 


Pu 2 (^ 2 ) = 


00 

J /Wafal) 3 ^) dx 1 . 


(2.7b) 



— oo 



Expected value (mean value): The quantity 

OO 

E{p) = J x Pn(x) dx (2.8) 

— OO 

is called the expected value (or mean value or statistical average) of the random variable 
/j, where E{-} denotes the expected value operator. The expected value operator E{-} is 
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linear, i.e., the relations E{ap} = aE{p} ( a £ IR) and E{p i + p 2 } = E{pi} + E{p 2 } 
hold. Let f{p) be a function of the random variable p. Then, the expected value of 
/(/i) can be determined by applying the fundamental relationship 



OO 

E{f(p)} = J f(x) P»(x) dx. 

— OO 



(2.9) 



The generalization to two random variables p\ and p 2 leads to 

OO OO 

E{f(pi,P 2 )}= J J f{xi,x 2 ) Pu 1 h 2 (xi,x 2 ) dxi dx 2 . 



— OO — OO 



(2.10) 



Variance: The value 



Var {p} = E{(n-E{^}) 2 } 

= E{n 2 }-{E{„}f (2.11) 

is called the variance of the random variable /z, where Var { • } denotes the variance 
operator. The variance of a random variable p is a measure of the concentration of p 
near its expected value. 

Covariance: The covariance of two random variables p\ and p 2 is defined by 



Cov {pi,p 2 } = E{(pi — E{pi}){p 2 — E{p 2 })} (2.12a) 

= E{pip 2 } — E{pi} ■ E{p 2 } . (2.12b) 

Moments: The kth moment of the random variable p is defined by 

OO 

E{p k } = J x k p^x) dx , k = 0,1,... (2-13) 

— OO 

Characteristic function: The characteristic function of a random variable p is 
defined as the expected value 

OO 

= E {e j2 ^} = J p^x) e j2 ™ x dx , (2.14) 

— OO 

where v is a real-valued variable. It should be noted that 'L M (— v) is the Fourier 
transform of the probability density function p^{x). The characteristic function often 
provides a simple technique for determining the probability density function of a sum 
of statistically independent random variables. 

Chebyshev inequality: Let p be an arbitrary random variable with a finite expected 
value and a finite variance. Then, the Chebyshev inequality 



(2.15) 




RANDOM VARIABLES 



15 



holds for any e > 0. The Chebyshev inequality is often used to obtain bounds on the 
probability of finding p outside of the interval E{p} ± ey/ Var {p}. 

Central limit theorem: Let p n (n = 1,2,..., iV) be statistically independent 
random variables with E{p n } = m (lii and Var {fin} = cr‘j ln . Then, the random variable 



A* = 



lim 

N — »oo 



l 

y/N 



N 

n— 1 



«V„) 



(2.16) 



is asymptotically normally distributed with the expected value E{p] = 0 and the 
variance Var {p} = a 2 = Jin^ ^ E^=i 

The central limit theorem plays a fundamental role in statistical asymptotic theory. 
The density of the sum (2.16) of merely seven statistically independent random 
variables with almost identical variance often results in a good approximation of the 
normal distribution. 



2.1.1 Important Probability Density Functions 

In the following, a summary of some important probability density functions often used 
in connection with channel modelling will be presented. The corresponding statistical 
properties such as the expected value and the variance will be dealt with as well. At 
the end of this section, we will briefly present some rules of calculation, which are of 
importance to the addition, multiplication, and transformation of random variables. 

Uniform distribution: Let 9 be a real-valued random variable with the probability 
density function 



Pe(x) 



1 r x 

— , X £ [— 7T, 7T) , 

0 , else . 



(2.17) 



Then, pe{x) is called the uniform distribution and 8 is said to be uniformly distributed 
in the interval [— n, 7 r). The expected value and the variance of a uniformly distributed 
random variable 9 are E{9} = 0 and Var {9} = 7t 2 /3, respectively. 

Gaussian distribution (normal distribution): Let p be a real- valued random 
variable with the probability density function 

(x — m ) 2 

p^x) = — =L — e 2<T " , x G 1R . (2.18) 

V27 roy, 

Then, p^x) is called the Gaussian distribution (or normal distribution) and p is said to 
be Gaussian distributed (or normally distributed). In the equation above, the quantity 
?n M £ IR denotes the expected value and er 2 £ (0, oo) is the variance of p, i.e. , 



E{p} = 



(2.19a) 
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and 



Var {/x} = E{h 2 } -m\ = o] l . (2.19b) 

To describe the distribution properties of Gaussian distributed random variables /x, we 
often use the short notation /i ~ 7V(?n M , a^) instead of giving the complete expression 
(2.18). Especially, for m ;t = 0 and a 2 = 1, N(0, 1) is called the standard normal 
distribution. 

Multivariate Gaussian distribution: Let us consider n real-valued Gaussian 
distributed random variables Hit^i ■ ■ ■ , Hn with the expected values m fli ( i = 
1,2, ...,n) and the variances (i = 1,2, ...,n). The multivariate Gaussian 
distribution (or multivariate normal distribution) of the Gaussian random variables 
Hi, H2, ■ ■ ■ , Hn is defined by 



PniU 2 .. Un ("Dj *^2) • • • j %n) — 



(V2tt) \/mc; 



o _i (x-mJ T C 



(2.20) 



where T denotes the transpose of a vector (or a matrix). In the above expression, x 
and are column vectors, which are given by 



( Xl \ 

X2 

elR” xl 

\ x n ) 



and 



( E{hi} \ 




( m Ml \ 


E{H2} 


= 


m V2 


V E{Hn} ) 




V / 



(2.21a) 



(2.21b) 



respectively, and det (C^ 1 ) denotes the determinant (inverse) of the covariance 
matrix 



C/i2A*l 


G/ 21 M 2 


•• C^ n \ 

CfJ. 2 Un 


e ET xra . 


(2.22) 




Cm 


■■ 







The elements of the covariance matrix C M are given by 

Cov {m, Hj} = E {(Pi ~ "G.X/H - m m)} - Vi, j = 1,2, (2.23) 

If the n random variables Hi are normally distributed and uncorrelated in pairs, then 
the covariance matrix C ^ results in a diagonal matrix with diagonal entries . . In 
this case, the joint probability density function (2.20) decomposes into a product of n 
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Gaussian distributions of the normally distributed random variables /z* ~ a 2 .). 

This implies that the random variables p $ are statistically independent for all i = 

1,2 



Rayleigh distribution: Let us consider two zero-mean statistically independent 
normally distributed random variables p \ and /z 2 , each having a variance (Tq, i.e. , 
p i,/i 2 ~ -ZV(0,<Tq). Furthermore, let us derive a new random variable from p i and /z 2 
according to £ = ^/zf + /z|. Then, £ represents a Rayleigh distributed random variable. 
The probability density function p^(x) of Rayleigh distributed random variables £ is 
given by 



Pc( x ) 




x > 0 , 

x < 0 . 



(2.24) 



Rayleigh distributed random variables £ have the expected value 



E{(} = vo^ (2.25a) 

and the variance 

Var{C} = a 0 2 (2~) . (2.25b) 



Rice distribution: Let /zi,/z 2 ~ N(0,ctq) and p £ IR. Then, the random variable 
£ = \J (/zi + p ) 2 + p 2 is a so-called Rice distributed random variable. The probability 
density function p^(x) of Rice distributed random variables £ is 



Pt, 0) 




x > 0 , 
x < 0 , 



(2.26) 



where Iq(-) denotes the modified Bessel function of Oth order. For p = 0, the Rice 
distribution p$(x) results in the Rayleigh distribution p^(x) described above. The first 
and second moment of Rice distributed random variables £ are [Wol83a] 



E{£} = cr 0 





and 




(2.27a) 



£{£ 2 } = 2a 2 + p 2 , 



(2.27b) 



respectively, where /„(•) denotes the modified Bessel function of nth order. From 
(2.27a), (2.27b), and by using (2.11), the variance of Rice distributed random variables 
£ can easily be calculated. 



Lognormal distribution: Let p be a Gaussian distributed random variable with 
the expected value and the variance cr 2 , i.e., /z ~ ./'/(m^, cr 2 ). Then, the random 
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variable A = e p is said to be lognormally distributed. The probability density function 
p\{x) of lognormally distributed random variables A is given by 



— — — P 2o V 7* f) 

VtoW ’ - ’ (2-28) 

0 , x < 0 . 

The expected value and the variance of lognormally distributed random variables A 
are given by 



E{ A} = 



(2.29a) 



and 

Var {A} = e 2m i+< (e< - l) , (2.29b) 

respectively. 



Suzuki distribution: Consider a Rayleigh distributed random variable ( with the 
probability density function P((x), according to (2.24), and a lognormally distributed 
random variable A with the probability density function p\(x), according to (2.28). Let 
us assume that £ and A are statistically independent. Furthermore, let 77 be a random 
variable defined by the product 77 = £ • A. Then, the probability density function p v (z) 
of 77 , that is 



Pv( z ) 



Z /'I 

\/27 rcr^CTu J y 3 

0 

0 , 



e 



2 « 2 °o 



_ (In a-m M ) 2 

e dy , z > 0 , 

z < 0 , 



(2.30) 



is called the Suzuki distribution [Suz77]. Suzuki distributed random variables rj have 
the expected value 

E{ri) = a QS J\e m ^% (2.31) 

and the variance 

Var {77 } = (Tq • e 2m '* +<7 M . ( 2e< - . (2.32) 



Nakagami distribution: Consider a random variable ui distributed according to the 
probability density function 



Pu{x) 



2 m m x 2m ~ 1 e~ ( m /^) x2 

0 , 



to > 1/2 , 
x < 0 . 



x > 0 , 



(2.33) 



Then, u denotes a Nakagami distributed random variable and the corresponding 
probability density function p^ (x) is called the Nakagami distribution or to- distribution 
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[Nak60] . In (2.33), the symbol T(-) represents the Gamma function, the second moment 
of the random variable lo has been introduced by Q = E{lo 2 }, and the parameter 
m denotes the reciprocal value of the variance of lo 2 normalized to O 2 , i.e., m = 
O 2 /E{(lo 2 — fl) 2 }. From the Nakagami distribution, we obtain the one-sided Gaussian 
distribution and the Rayleigh distribution as special cases if m = 1/2 and m = 1 , 
respectively. In certain limits, the Nakagami distribution, moreover, approximates 
both the Rice distribution and the lognormal distribution [Nak60, Clra79]. 

2.1.2 Functions of Random Variables 

In some parts of this book, we will deal with functions of two and more random 
variables. In particular, we will often make use of fundamental rules in connection with 
the addition, multiplication, and transformation of random variables. In the sequel, 
the mathematical principles necessary for this will briefly be reviewed. 

Addition of two random variables: Let /j-| and be two random variables, which 
are statistically characterized by the joint probability density function Pmm( x 1 , 2 : 2 ). 
Then, the probability density function of the sum p = /ii + /i 2 can be obtained as 
follows 

OO 

Pn(y) = J - x-i) dx x 

— OO 
OO 

= J Pmn 2 (y - X 2 ,X 2 ) dx 2 ■ (2.34) 

— OO 

If the two random variables pi and P 2 are statistically independent, then it follows 
that the probability density function of p is given by the convolution of the probability 
densities of pi and p 2 - Thus, 



p»(y) = Pm(y)*p^(y) 

00 

= J PnA x i)p^{y-x\) dx! 

— OO 
OO 

= J Pm {y - X2)Pu , 2 (x 2 ) dx 2 , (2.35) 

— OO 

where * denotes the convolution operator. 

Multiplication of two random variables: Let ( and A be two random variables, 
which are statistically described by the joint probability density function p^\(x, y). 
Then, the probability density function of the random variable 77 = £ • A is equal to 



OO 




— OO 



(2.36) 
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From this relation, we obtain the expression 




■P\{y) dy 



(2.37) 



for statistically independent random variables £, A. 



Functions of random variables: Let us assume that H\, fi 2 , ■ ■ ■ , jx n are random 
variables, which are statistically described by the joint probability density function 
Pmn^i.iJ. n ( x i, x 2 , • • • , £«)• Furthermore, let us assume that the functions /i, fi, ■ ■ ■ , fn 
are given. If the system of equations fi(x 1 , £ 2 , . . . , x n ) = Vi (i = 1,2,..., n) has real- 
valued solutions £i„, X 2 v, ■ ■ ■ , x nv (u = 1,2, . . , , m), then the joint probability density 
function of the random variables £1 = fi(pi, P 2 , ■ ■ ■ , Pn), £2 = f 2 ( 1 x 1 , 1 x 2 , ■■■ , fx n ), ■ ■ ■ , 
£n = fn(/x 1 , 1 x 2 , ■■■, (Xn) can be expressed by 



m 

Pii£2...tn(yi>y2,--->yn) = 

V=1 

where 



Pm 



X x i vi %2vi • • • 1 %nv) 



| Ji^lvi %2 vi • • • 1 %nv) | 



J(x l,x 2 , ...,x n ) = 



M 


dfi . 


.. Ah 


dx\ 


dX2 


dx n 


dh 


dh . 


.. 


dxi 


dx 2 


dx n 



dh 


df n . 


. . ih. 


dx\ 


6x2 


dx n 



denotes the Jacobian determinant. 



(2.38) 



(2.39) 



Furthermore, we can compute the joint probability density function of the random 
variables £ 1 , £ 2 , . . . , for k < n by using (2.38) as follows 



Pii&.jk (vi > y2 , 




00 

■ J Ptit2...£ n (yi>y2, ■ ■ ■ ,y n ) dy k +i dy k+2 ...dy n . 

— OO 

(2.40) 



2.2 STOCHASTIC PROCESSES 

Let (Q,A,P) be a probability space. Now let us assign to every particular outcome 
s = Sj £ Q of a random experiment a particular function of time fi(t, si) according to 
a rule. Hence, for a particular .s, £ Q, the function y(t, si) denotes a mapping from IR 
to IR. (or C) according to 



I x(-,Si ) : IR — > IR (or C) , 1 1 — > y(t, si) . 



(2.41) 
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The individual functions p(t, si) of time are called realizations or sample functions. A 
stochastic process /.i(t, s ) is a family (or an ensemble) of sample functions p(t, sf), i.e. , 
p(t,s) = {p,(t,Si)\si G Q} = 

On the other hand, at a particular time instant t = to G IR, the stochastic process 
fi(to,s) only depends on the outcome s and, thus, equals a random variable. Hence, 
for a particular to G IR, p(to, s) denotes a mapping from Q to IR (or C) according to 

p(t 0 , •) : Q -> IR. (or C) , s /t(f 0 , s) . (2.42) 

The probability density function of the random variable /t(fo, s) is determined by the 
occurrence of the outcomes. 

Therefore, a stochastic process is a function of two variables t G IR and s G Q, so that 
the correct notation is /i(t, s). Henceforth, however, we will drop the second argument 
and simply write /i(t) as in common practice. 

From the statements above, we can conclude that a stochastic process /x(t) can be 
interpreted as follows [Pap91]: 

(i) If t is a variable and s is a random variable, then /i(t) represents a family or an 
ensemble of sample functions p(t, s). 

(ii) If t is a variable and s = So is a constant, then p(t) = /i(t, Sq) is a realization or 
a sample function of the stochastic process. 

(iii) If t = to is a constant and s is a random variable, then n(to) is a random variable 
as well. 

(iv) If both t = to and s = so are constants, then p(to) is a real-valued (complex- 
valued) number. 

The relationships following from the statements (i)-(iv) made above are illustrated in 
Figure 2.1. 



Complex-valued stochastic processes: Let //(f) and //'(f) be two real-valued 
stochastic processes, then a (complex-valued) stochastic process is defined by p(t) = 
//(t) + 

We have stated above that a stochastic process /t(f) can be interpreted as a random 
variable for fixed values of t G IR. This random variable can again be described by 
a distribution function /’^(x; f) = P(/t(f) < x) or a probability density function 
p^,(x;t) = dF^(x',t)/dx. The extension of the concept of the expected value, which 
was introduced for random variables, to stochastic processes leads to the expected 
value function 

m^(t) = E{p(t)} . (2.43) 

Let us consider the random variables p{t\) and p,(t 2 ), which are assigned to the 
stochastic process p{t) at the time instants t\ and t 2 , then 



r^(ti,t 2 ) = E{p?( fi)/i(f 2 )} 



(2.44) 




